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ABSTRACT 

The  terms  of  Sutcliffe»s  development  formula  are  evalu- 
ated by  finite  differences  for  one  hundred  cyclones.   The 
values  of  relative  divergence  obtained  are  correlated  with 
the  subsequent  24-hour  pressure  change  and  24-hour  average 
speed,   A  correlation  of  -.08  was  obtained  between  the  relative 
divergence  at  the  center  and  the  subsequent  24-hour  pressure 
change.   Computing  the  algebraic  sum  of  relative  divergence 
across  the  center  and  correlating  this  value  with  the  24-hour 
pressure  change  gave  a  coefficient  of  .28.   Using  the  alge- 
braic difference  of  relative  divergence  across  the  center 
and  correlating  this  with  the  subsequent  24-hour  average 
speed  gave  a  coefficient  of  -.02.   It  therefore  appears  that 
the  finite  difference  evaluation  of  the  Sutoliffe^  develop- 
ment equation,  as  made  in  this  investigation,  is  not  a  useful 
method  for  a  24-hour  prognosis  of  the  deepening,  filling,  or 
movement  of  winter  cyclones.  , 
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1.   Introduction. 

Any  attempt  to  forecast  the  future  position  or  intensity 
of  a  cyclonic  center  should  be  based  on  the  present  dynamic 
processes  as  well  as  the  changes  which  these  processes  will 
undergo.   The  latter  part  of  this  statement  involves  the  more 
difficult  phase  of  forecasting.   There  have  been  various 
approaches  to  the  solution  of  this  problem,  none  have  been 
one  hundred  percent  effeotive.   In  some  instances  the  fore- 
oast  is  entirely  subjective  and  depends  primarily  on  the 
experience  of  the  forecaster.   On  the  other  hand,  empirical 
or  statistical  rules  have  been  devised  which  give  satisfactory 
results  in  some  geographical  areas.   It  was  the  purpose  of 
this  Investigation  to  see  if  Sutcliffe's  development  formula 
r3"j  would  serve  as  an  accurate  indication  of  the  future  in- 
tensity of  a  cyclonic  center  and  the  average  speed  it  would 
maintain  for  a  twenty-four  hour  period. 


2*   Resume  of  the  derivation  of  Sutcliffefs  development 
formula* 

It  may  be  desirable  to  briefly  outline  the  derivation 

of  this  development  formula  and  note  the  assumptions  that 

are  made*  Essentially  he  begins  with  the  general  vorticity 

equation  in  the  form 

where  pressure  is  taken  as  the  vertical  variable*  According 
to  estimates  by  Sutoliffe,  the  last  term,  called  the  "tilt" 
term,  is  largely  offset  by  ^(^p'^/)  •     Another  modification 
is  made  by  neglecting  the  vertical  advection  of  vorticity 
which  is  normally  believed  to  be  small*   These  simplifications 
reduce  equation  (1)  to 

Equation  (2)  involves  the  divergence  at  any  one  level,  but 
the  local  field  of  pressure  is  not  the  result  of  the  motion 
at  one  level,  but  rather  of  the  integrated  effects  of  the 
three  dimensional  motion*   In  order  to  get  an  estimate  of 
the  contribution  of  a  layer,  the  difference  is  taken,  result- 
ing in 

f (v, ■  K  -  vf  •  v,.;  -  UkD  +  u.-^.*rt  -  v*  •  v,tf <  f)   (si 

«»  v. 

where  the  subscript  "0"  denotes  the  lower  layer*  Applying 
the  geostrophic  approximation  for  vorticity  in  the  form  of 
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r-f  v;* 


the  local  rate  of  change  of  the  difference  in  vorticity 
becomes 

where  "h"  is  the  thickness  between  the  two  levels.   Substi- 
tution of  equation  (4)  into  equation  (3)  gives 

f  (.v,  •  v,  -  v,-v„8 )  =  -  j  v;  g  * V„0  •  yP  ($> f )  -  v„  •  Vf ♦  H .  (5) 

In  order  to  use  this  equation  as  a  prognostic  tool,  the 
term  of  4^?  must  be  estimated  or  replaced  by  some  reasonable 
assumption  for  the  thickness  change.  .  In  this  derivation  it 
is   assumed  that  temperature  changes  and  therefore  thickness 
changes  are  purely  advective;  and,  finally,  that  this  advection 
may  be  approximated  by  the  geostrophic  wind.   That  is  ^r-V^*\|h. 
With  this  assumption,  non-adiabatic  effects  and  vertical 
motions  are  neglected,  and  in  some  areas  these  may  have  marked 
affects  on  the  thermal  field.  Using  the  advective  assumption, 
as  well  as  the  equations,  Vj  =  V0  * V^    and  f    ~   Jj6t  fT  $ 

the  final  form  of  equation  (5)  gives  the  "relative  divergence" 
of  Sutcliffe  as 

(ol)      0>)  U> 

where  s1  is  the  arc  length  in  the  direction  of  the  thickness 

lines. 

We  will  now  examine  each  of  the  various  terms  of  this 

equation. 

a.  From  term  (a)  it  is  indicated  that  thermal  wind  over 
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a  surface  center  of  positive  vortioity  implies  relative 
divergence  ahead,  relative  convergence  to  the  rear.  However, 
this  term  has  no  effect  on  the  intensification  of  the  center, 
because  at  the  center  ~?   is  zero.   Hence  it  merely  affects 
the  propagation  of  the  wave. 

*>•  ""  r  "W    This  term  indicates  that  there  should  be 
relative .divergence  where  the  thermal  vorticlty  decreases 
in  the  direction  of  the  thermal  wind  and  relative  convergence 
where  the  thermal  vorticity  increases  In  the  direction  of 
wind  shear.   In  the  normal  synoptic  situation,  where  the 
thickness  trough  Is  In  the  rear  of  a  surface  depression  and 
a  thickness  ridge  over  or  ahead  of  it,  there  would  be  a  de- 
crease of  thermal  vorticity  downstream  over  the  surface  de- 
pression, which  would  mean  relative  divergence.   On  the  other 
hand,  when  the  thermal  trough  precedes  the  surface  depression, 
there  should  be  relative  convergence. 

c.  -^r  ltt  Prom  this  term  one  would  get  relative 

convergence  when  the  thermal  wind  has  a  poleward  component 
and  relative  divergence  with  an  equatorward  component.   How- 
ever, in  mid-latitudes,  the  magnitude  of  the  rate  of  change 
of  the  coriolis  parameter  is  generally  small  compared  to  the 
other  terms  and  for  purposes  of  this  Investigation  will  be 
neglected. 

It  has  been  observed  that  divergence  varies  vertically 
and  that  the  normal  distribution  over  a  surface  low  system 

is  to  have  convergence  in  the  lower  layers  with  divergence 
above.   In  order  to  have  development,  the  integrated  divergence 


aloft  must  over  compensate  for  the  integrated  convergence  in 
the  lower  layers.   The  amount  of  this  over  compensation 
should  be  a  measure  of  the  deepening  or  filling  of  the  sur- 
face system.   Instead  of  measuring  the  divergence  and  con- 
vergence for  the  entire  vertical  extent  over  a  system, 
measuring  it  at  some  specified  level  or  levels  might  be 
feasible  and  still  give  qualitatively  similar  results. 
Assuming  a  schematic  vertical  distribution  of  divergence 
over  a  cyclone  or  anticyclone  as  in  Figure  1,  Sutcliffe 
implies  that  the  differences  of  divergence  of  any  two  levels 
should  be  an  indication  of  the  deepening  or  filling  of  the 
surface  system. 
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Figure  1 

Schematic  representation  of  the  distribution  of  divergence. 


That  is,  if  the  relative  divergence  over  a  cyclone  is  a  large 
positive  value,  the  deepening  should  be  large  and  as  this 
value  decreases,  the  deepening  should  also  decrease.  An 
equivalent  statement  for  anticyclones  would  associate  negative 
values  of  relative  divergence  with  anticyclonic  development. 


3#   Investigation  procedures  and  correlations. 

The  purpose  of  this  investigation  was  to  see  if  it ' 
would  be  possible  to  correlate  the  present  values  of  relative 
divergence  as  measured  by  terms  (a)  and  (b)  of  equation  (6) 
with  a  twenty-four  hour  deepening  or  filling  of  a  surface 
low,  as  well  as  the  average  speed  of  movement  over  the  same 
period  of  time.   One  hundred  cases  of  surface  lows  in  the 
period  of  December  1956  through  February  1957  were  used  for 
the  investigation.   These  lows  were  located  in  mid-latitudes, 
ranging  from  28°N  to  60°N.   They  were  selected  primarily  over 
the  North  American  continent  in  order  to  have  the  advantage 
of  a  dense  reporting  network*  thereby  reducing  the  analysis 
error  due  to  subjectivity. 

The  initial  attempt  in  the  investigation  was  to  evaluate 
the  relative  divergence  at  the  center  of  a  surface  low.   In 
regard  to  Sutoliffe's  formula  this  meant  evaluating  the  term 

-XL    ilr     .   The  other  term  would  have  no  contribution  because 
the  surface  vorticity  center  and  low  center  are  normally  co- 
incident, thereby  making  Ale,   equal  to  zero.   Positions 
three  hundred  nautical  miles  upstream  and  downstream  along 
the  same  thickness  line  were  used  and  the  calculation  of 
thermal  vorticity  at  these  points  was  made  by  the  finite 
difference  approximation  of   ^2  \7  h  $   namely, 

±2  ( K  -  h \\  .   a  grid  distance,  d,  of  5  degrees  latitude 

was  used  which  is  consistent  with  the  values  used  in  graphical 
and  numerical  techniques.   A  second  finite  difference  approx- 
imation was  used  to  evaluate  the  rate  of  change  of  thermal 


vortioity  across  the  center.   The  latitude  of  the  surface 
low  center  and  the  geostrophic  thermal  wind  over  the  center 
were  the  other  values  used  in  the  evaluation, 

A  plot  of  the  relative  divergence,  arrived  at  by  this 
method,  versus  the  twenty-four  hour  pressure  change  is  shown 
in  Figure  2.   There  is  no  definite  grouping  of  the  values, 
which  indicates  that  the  future  deepening  or  filling  of  the 
system  cannot  be  predicted  using  the  relative  divergence  at 
the  center.   A  correlation  coefficient  of  -.08  was  computed 
from  the  data  of  this  evaluation. 

The  results  of  the  initial  investigation  were  not  con- 
clusive so  another  approach  was  taken.   In  accordance  with 
the  kinematic  formula  of  deepening  or  filling  of  the  surface 
center  should  be  proportional  to  the  algebraic  sum  of  the 
divergence  ahead  and  behind  the  system  ClJ  •   That  is,  if 
there  is  more  relative  divergence  ahead  of  the  system  than 
there  is  relative  convergence  behind  the  system,  the  system 
should  deepen.   To  test  this,  it  was  necessary  to  evaluate 
terms  (a)  and  (b)  of  Sutcliffe^  development  formula.   The 
values  of  thermal  vorticity  and  surface  vorticity  for  positions 
along  the  same  thickness  line  six  hundred  nautical  miles  ahead 
and  behind  the  surface  low  center,  as  well  as  at  the  center 
were  calculated.   The  geostrophic  finite  difference  approxi- 
mations of    -^  y7  h  and    £  ^J  &a  were  used. 
The  finite  difference  approximation  was  used  to  compute  the 
rate  of  change  of  these  quantities  along  the  thickness  lines 
for  points  three  hundred  nautical  miles  ahead  and  behind  the 
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center.   The  thermal  wind  and  latitude  of  the  points  three 
hundred  nautical  miles  ahead  and  behind  the  surface  center 
were  the  other  values  used  in  the  computation.  The  algebraic 
sum  of  terms  (a)  and  (b)  were  computed  across  the  low  and 
plotted  versus  the  actual  deepening  or  filling  for  twenty- 
four  hours*  Figure  3  shows  this  plot  and  the  scatter  of 
this  diagram  implies  that  very  little  can  be  inferred  from 
the  present  values  of  relative  divergence  sums  and  the  subse- 
quent pressure  change,  A  correlation  coefficient  of  #28 
between  the  relative  divergence  sums  and  the  twenty-four 
hour  pressure  change  was  computed  from  the  data  in  this  test. 

The  third  phase  of  the  investigation  consisted  of  an 
attempt  to  correlate  the  relative  divergence  with  the  twenty- 
four  hour  speed  of  the  surface  low  system.  The  speed  of  the 
system  should  be  proportional  to  the  gradient  of  relative 
divergence  across  the  system.   To  evaluate  this  quantity, 
the  algebraic  difference  of  the  relative  divergence  values 
for  the  points  three  hundred  nautical  miles  ahead  and  behind 
was  taken.  The  subsequent  speed  of  the  associated  low  was 
evaluated  in  units  of  degrees  of  latitude  per  twenty-four 
hours,  A  scatter  diagram  of  the  speed  versus  the  difference 
of  relative  divergence  is  shown  in  Figure  4,   The  correlation 
coefficient  of  these  two  quantities  is  -,02,  and  thus  no 
definite  statement  can  be  made  as  to  what  average  speed  can 
be  expected  from  the  value  of  the  difference  of  relative 
divergence  across  the  system. 
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4.   Results  and  interpretations. 

As  a  result  of  the  observations  and  computations  made 
in  this  investigation,  it  appears  that  on  the  basis  of 
Sutoliffe»s  development  formula  the  future  development  or 
displacement  cannot  be  forecast. 

The  poor  relationships  resulting  from  the  various  phases 
of  the  investigation  could  be  attributed  to  the  fact  that 
the  actual  vertical  distribution  of  divergence,  as  well  as 
relative  divergence  is  more  complex  than  that  which  was 
assumed  by  Sutcliffe.  Also  some  of  the  assumptions  in- 
volved in  the  derivation  of  the  development  formula  are  not 
completely  realized  at  all  times.  Perhaps  the  length  of 
time  that  was  chosen  for  the  future  Intensity  was  too  long, 
but  on  a  cursory  examination  of  the  twelve-hour  pressure 
changes  and  speed,  no  Improvement  was  apparent.  It  also 
may  be  indicated  that  the  distances  used  in  the  evaluation 
may  be  too  large  and  that  positions  closer  to  the  center 
would  give  better  approximations  for  the  finite  difference 
method.   Perhaps  making  some  modifications  in  the  method  of 
evaluation,  or  by  empirically  allowing  for  the  neglected 
terms  of  the  development  formula,  a  more  promising  result 
could  be  obtained.  An  evaluation  of  the  terms  of  Sutcliffe !s 
development  formula  on  a  synoptic  chart  by  use  of  a  small 
grid  system  and  Taylor  Series  expansion  has  been  carried  out 
by  Sawyer  and  Matthewman  (_2j/ ',  but  no  attempt  was  indicated 

of  correlating  the  results  with  the  future  intensity  of  the 
cyclonic  system. 
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